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Intriguing quantum effects that result from entangled molecular rovibrational states are shown to 
provide a novel means for controlling both differential and total collision cross sections in identical 
particle diatom-diatom scattering. Computational results on elastic and inelastic scattering of para- 
H2 and para-H2 are presented, with the collision energy ranging from 400 cm -1 to the ultracold 
regime. The experimental realization and possible extension to other systems are discussed. 



I. INTRODUCTION 
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^ \ Recent years have witnessed an increasing interest in controlling molecular and atomic processes jl| The 
/"N ■ central principle of Coherent Control (CC) teaches that such control can be achieved by using coherent laser fields to 
induce constructive or destructive interference between multiple indistinguishable pathways to the same final state. 
The yield of desired or undesired product can thus be enhanced or suppressed by manipulating the resultant quantum 
interference terms. Thus far, most experimental and theoretical studies of CC of atomic and molecular processes have 
focused on unimolecular reactions or half collisions, such as photodissociation, photoionization and photoassociation. 
By contrast, controlling full collision processes such as atom-atom scattering, atom-molecule scattering, and 
, molecule-molecule scattering remains a significant challenge. Existent studies along this direction have largely em- 
ployed a number of approaches [Q, [|, [| |?], || . The first ^, ||, |(| actively manipulates coherent light sources to alter 
the interaction potential of the colliding particles. For example, in ultracold atom-atom collisions, the scattering cross 
sections can be extensively modified by using laser fields that are in near resonance with one of the excited molecular 
electronic states [Q, or by employing a d.c. field that induces an additional dipole-dipole interaction potential 
However, it appears difficult to use this approach to control atom-molecule or molecule-molecule scattering, since the 
number of degrees of freedom is much larger than that in atom-atom scattering. A second approach is laser catalysis 
H in which the laser field photons serve as a catalyst. In this case the laser accelerates (or suppresses) the reaction 
by inducing a barrier hopping process, via a virtual excitation to a bound excited state, while remaining unchanged 
at the end of the event. The third approach, which is the subject of this study, is coherent control j7]. Although the 
general theory of collisional CC has been worked out [Q and conditions for control established || the experimental 
problem remains as to how one can introduce the required adjustable quantum phases and amplitudes into coherent 
scattering processes in a fashion so as to introduce control based upon quantum interference. That is, several recent 
studies on CC of collisions in model scattering systems S make clear that, for collisions of non-identical molecules, 
coherent control requires that one establish subtle quantum correlations between the translational and internal degrees 
of freedom of the reactants. Setting up the required initial states presents an experimental challenge, since it requires 



coherent matter-wave beams in which the internal states of the molecule are correlated in a precise way with the 
- - 1 translational states. 

Such correlated translational-internal states are not required, however, if one attempts coherent control in the 
collision of identical molecules, i.e., AB + AB. Rather, it suffices [nil to prepare the reactants in superpositions of 



initial rovibrational states. In this case, in accord with Ref. |10 , the collision cross sections can be controlled by 
manipulating the character of the initial coherent superposition states in which the molecules are prepared. 

In this paper we examine the coherent control of identical particle scattering in detail, with para-H2 + para-H2 
as an example. In doing so we approach this issue from a different perspective from that in our earlier work p0[ , 
demonstrating deep connections to the construction and use of entangled states, a topic of great recent interest in 
quantum information science That is, we show that the scattering amplitudes, and the resultant differential and 
total collision cross sections, depend strongly on quantum entanglement embedded in the initial rovibrational state. 
In essence, entanglement-induced quantum effects are shown to occur without the explicit preparation of entangled 
states. This provides a novel means of controlling both the the differential and total collision cross sections in identical 
diatom-diatom scattering. 



* Permanent address: Chemical Physics Department, The Weizmann Institute of Science, Rehovot, Israel 76100 



2 



This paper is organized as follows. In Sec. 0, we introduce some concepts related to quantum entanglement 
between identical particles and then address the theory of identical diatom-diatom scattering with particular reference 
to entangled molecular rovibrational states as the incoming asymptotic states. The theoretical formalism is presented 



Sec. Ill for para-H2 + para-H2 system, without loss of much generality. Results for the scattering of this system 
are shown in Sec. |^ where detailed computational results on the role of entangled rovibrational states in both elastic 
and inelastic scattering, with the collision energy ranging from 400 cm -1 to the ultracold regime, are explored. A 
discussion is provided in Sec. 



II. AB + AB SCATTERING 



A. Identical Particle Entanglement 



Quantum entanglement is a striking feature of quantum physics. Qualitatively, two particles are said to be quantum 
entangled if their total wavefunction is inseparable. In particular, if two particles are entangled then neither of them 
possesses a complete set of properties, since measuring one of them would collapse the total wavefunction and thus 
affect the properties of the other. Quantum entanglement is a fundamental issue of considerable theoretical interest, 
and plays a key role in various modern areas such as quantum teleportation ]l2| , quantum cryptography, and 
quantum computing plj| . 

In the case of identical particles the total wavefunction, including all degrees of freedom, is always symmetrized or 
anti-symmetrized with respect to particle exchange. Hence, even if the total wavefunction of two identical particles is 
simply obtained by symmetrizing (or anti-symmetrizing) a separable state, it appears to be unfactorizable, and hence 
apparently entangled. This is, in fact, misleading and hints that quantum entanglement between identical particles is 
somewhat more subtle than quantum entanglement between non-identical particles. In a recent discussion, Ghirardi 
et al [fhlf have given clear definitions of quantum entanglement between two identical particles: (i) two identical 
fermions are entangled if the total wavefunction cannot be obtained by anti-symmetrizing a factorized state, and (ii) 
two identical bosons are entangled if the total wavefunction cannot be obtained by symmetrizing a factorized product 
of two orthogonal states and if the total wavefunction is not a product of the same state for the two particles. Under 
these circumstances, the state is entangled and each of the particle pairs does not possess a full set of properties. 

One intriguing aspect of quantum entanglement between identical particles, relevant to the discussion below, is that 
it can force different permutation symmetries on some degrees of freedom. For example, consider a system of two 
freely moving identical spin 1/2 particles, whose spin up and spin down states are represented by | "fi), | ji), | T 2 ), 
I | 2 ). If the total spin is zero, then the two spins are entangled and the resulting ("EPR state") \ip) is given by 

M = ^(|ti>®Ua)-Ui>®|T2». (1) 

Note that \tp) in Eq. (Q) cannot be obtained by anti-symmetrizing a factorizable total wavefunction (since such an 
anti-symmetrized state should also describe the translational motion) of two spin 1/2 particles. Clearly, in Eq. 
([j]) acquires a factor of — 1 upon permutation. Since the permutation symmetry for the total wavefunction including 
all degrees of freedom is —1 for identical fermions, the permutation symmetry of the spatial degrees of freedom must 
be +1. Thus, in a simple bound system such as para-H2, the rotational quantum number has to be even. Likewise, in 
unbound cases such as electron-electron or proton-proton scattering, the differential cross sections do depend strongly 
on how the two spins are entangled, even if the scattering potential is spin independent, e.g., a pure Coulomb potential 
|^5| . This selection of specific permutation symmetry combinations due to entanglement of a subset of the system 
degrees of freedom is central to the results below. 



B. The Role of Entangled Rovibrational States in Identical Particle Scattering 



Consider identical particle scattering. We focus on the diatom-diatom case AB + AB, but the considerations are 
general. If | j, m, v) denotes an eigenstate of the AB diatom with angular momentum quantum number j, angular 
momentum projection quantum number m and vibrational quantum number v, then a typical entangled ro- vibrational 
state of AB + AB is of the form 



{r 1 ,r2\ip(jim 1 v 1 j2m 2 V2)) a ,p = cos(a)(ri|ji, mi, Ui)(r 2 |j 2 , m 2 , v 2 ) 

+ sin(a) exp(i/3) (r 1 \j 2 , m 2 , v 2 ) (r 2 mi , vi) . 



(2) 
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Of specific interest later below are the entangled states 

\ip(jimiV 1 j 2 m 2 V2))± = -^={(r 1 \j 1 ,m 1 ,v 1 }(r 2 \j2,m2,V2} 

±(ri\j 2 ,m2,V2)(r 2 \j 1 ,m 1 ,v 1 )), (3) 

where the only distinction between the two states |V , 0'i m i' l, i.72 m 2^2))± is the relative phase (3 (0 or 7r) between the 
two participating states (Yi\jiniivi) {Y 2 \j 2 m2V2) and (ri|j2m2V2)(r2|jiTOiVi). This difference is significant insofar 
as \ip(jimiVij 2 m 2 V2)} + is invariant upon permutation of ri and r2, whereas \ip(jimiVij 2 m2V2))- acquires a fac- 
tor of —1 upon permutation of ri and T2- ft should be stressed that the permutation symmetry associated with 
\4'(ji Tn i v ij2m2V2))± has nothing to with the permutation symmetry of the total wavefunction of identical particles, 
since the states |'0(ii TO i' l 'i.72W2i>2))± describe only rotational and vibrational motion. Note also, for use later below, 
that entangled states \ip(ji m i v ij2 m 2 v 2))± & lso form a new set of basis states to describe a general entangled state, 
i.e., 

. ^> cos(a) +sin(a)exp(i/3) . 
\ip(jim 1 vij2m2V2)) a ,f3 = \tp{JimiVij2m 2 V2j) + 

cos(a) -sin(a)exp(z/3) . 
H \^{jim 1 v 1 j 2 m 2 V2j)-- (4) 

Suppose that prior to the AB + AB collision the AB molecules (here labeled AB and AB' for convenience) are 
prepared, in accord with Ref. pi, in the superposition states: 

\ipAB) = cos(ai)|jiTOiwi) + sin(ai) exp(i/3i)|j 2 TO2f2) 

\ipAB') = cos(a 2 )|jiTOiwi) + sin(a 2 ) exp(i/3 2 )|j2m2W2) (5) 

Then the total internal wavefunction (before symmctrization) is a direct product \ipdp) of these two superposition 
states, rather than an entangled molecular state. Nevertheless, \ipdp) can be expressed in terms of an entangled state 
plus two additional components, i.e., 

\ipd P ) = yexp(i/3 2 )|^(jimiUij 2 TO2W2))a/3 + cos(ai) cos(a 2 )|jimiwi) <gs \jimiVi) 

+ sin(ai)sin(a 2 )exp[i(/3i + /3 2 )] \j2m2V2) ® |j 2 m 2 u 2 ), (6) 

where the entangled state component |V'0'i m i l, ij2"^2^2)) a/3 is given by Eq. (|]), with 



y = ycos 2 (ai)sin (a 2 )+sin (ai) cos 2 (a 2 ), 
a = cos _1 [cos(ai) sin(a 2 )/y], 

= Pi- 02, (7) 

and the two additional components \jimiVi) g> [i imivi) and \j2m2V2) ® ^21^2^2) are called, in the spirit of previ- 
ous coherent control work, "satellite" states [[[(J. From this viewpoint, the contribution from the entangled state 
component is the control term. 

If our interest is in coherent control of AB + AB scattering, then this provides an obvious route to control. That is, 
by altering the a.; and (3i in the prepared diatomic state, one can alter the initial state. This introduces controllable 
quantum interferences, as previously discussed [ [To| , hence altering product cross sections. 

Alternatively, if our interest is in the dynamics of entangled states then collisions of AB + AB in the prescribed 
initial states may well provide a route to directly observing entanglement-induced quantum effects. Specifically, there 
clearly are cases, some cited in computations below, where contributions from the entangled molecular state dominate 
over that of the satellite states. That is, for some channels characterized by j[, v[, j' 2 , v' 2 , the scattering amplitude 
for the initial direct product state \ipdp) is entirely due to the y e^jp(i/3 2 )\ip(jimiVij2m2V 2 )) a /3 component. For these 
cases, we have 

\ipdp)) ~ V 2 v (0, \^{jim 1 v 1 j2m2V 2 ) al3 )) , (8) 

where a(9, \ip)) represents the differential cross section summed over m! x and m 2 and integrated over the azimuthal 
angle <f> for the initial state \tp). Alternatively, and more difficult to implement, if we restrict the measurement to a 
particular total energy of the product, contributions of the satellite states can be zero due to the energy restriction 
and the scattering becomes essentially that due to an initially entangled molecular state | - 0(jimiUij2W2V2))a/3- Thus, 
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significantly, in either case one can experimentally observe and utilize entanglement-assisted coherent control effects 
in diatom-diatom scattering without the difficult task of preparing entangled molecular states. 

Below we develop a general description of the quantum effects associated with entanglement of molecular rovibra- 
tional motion in bimolecular scattering. We focus on the contribution from the entangled state since the remainder 
(i.e., the satellite terms) is ordinary scattering. In doing so we specifically study the scattering system para-H2 + 
para- H2, the simplest case of diatom-diatom scattering. Although many of the following arguments apply to both 
reactive and nonreactive scattering, we confine attention to elastic and inelastic scattering only. 



III. SCATTERING FROM ENTANGLED RO VIBRATIONAL STATES: FORMALISM 



Consider AB + AB scattering. We represent the vector pointing from the center of mass of the projectile to 
that of the target by R, and the relative inter-atomic distance vectors of the two para-H2 molecules by ri and r2 , 
respectively. We use ji,rm,Vi (i = 1,2) to denote the quantum numbers of angular momentum, the projection of 
angular momentum onto a space-fixed z axis, and vibrational motion, respectively, for each para-H 2 molecule. The 
fact that both the nuclear spin and the electronic spin of para-H2 molecules are zero greatly simplifies the problem, 
while retaining the essence of physics. Note that each para-H2 molecule (with electronic ground state) is a boson, so 
that the permutation symmetry of the total wavefunction |^(ri, r 2 , R)) in the scattering problem is +1. That is, 

|*(r 1 ,r 2 ,R)) = |*(r 2 ,r 1 ,-R)). (9) 

As we shall see more clearly below, this implies that the parity of |^(ri, ^2, R)) with regard to R depends on the 
permutation symmetry of the internal degrees of freedom ri and T2- 

In traditional scattering, the initial state is characterized by fc,z, jx,mi,vx, J2 3 rn -2 •> ^2 j and the final state charac- 
terized by £/,R, j[, m^v^, j' 2 ,m 2 ,v' 2 , where fik and tik 1 are the translational momenta for the initial and final states, 
with k = k! (k 7^ k') in the case of elastic (inelastic) scattering. Below we exclude the case of (ji, m\,vi) = (j2, m 2 , V2) 
as it does not allow for entanglement in rovibrational states. Since the colliding molecules are identical, both the 
incoming and the outgoing asymptotic states have to be appropriately symmetrized. However, to obtain the scattering 
amplitude it suffices to symmetrize either the incoming or outgoing state only |lq |. If we choose to symmetrize the 
incoming state then the scattering amplitude is given by 

f(k''Rj' 1 m' l v' l 3 2 m 2 v 2 <- kij 1 m 1 v 1 j 2 m2V2) 
= f(k'*k,j' 1 m' l v' 1 j'2m 2 v 2 <— kijimiVij 2 m 2 v 2 ) 

+.f(k'Ilj' 1 m' 1 v' 1 j 2 m 2 v 2 <- k(-z)j 2 m 2 v 2 jim 1 v 1 ), (10) 

where / represents the unsymmetrized scattering amplitudes under the assumption that the two molecules are dis- 
tinguishable. Specifically, in terms of the T-matrix elements, we have [|l6| [l7) 

/(fc'RjJmXij2 TO 2 w 2 *~ kijim 1 v 1 j2m 2 V2) 



7 E E E ^ w '^(R)« J>i;2 C mi2 



and 



Vkk 

JMm' lj\ 2 m\ 2 l'j' 12 m' 12 

X ^>i^^ m2 TJM (j^^^2j;2^bl^j2.2j 1 20, (11) 



f{k''R.j' 1 m' 1 v' 1 j 2 m'2v' 2 <- k(-z)j 2 m2V2jim 1 v 1 ) 



7 E E E (-l) ^ ^^m^ ^ >™'(R)^ ;2mi2 Q^ 2ml2 



JMm' lji2m 12 l'j' 12 m[ 

-li'v2. m 'v2 /Oj'l2"ll2 rpJM ( ■/ „,/ •/ „,/ „•/ // 



"(^^(^^"VM&mvihvd^ (12) 



where Y™ is the spherical function and Cjimij 2 m 2 are tne Clebsch-Gordan coefficient. 

Consider now the case where the molecules are initially prepared in an entangled state of molecular rovibrational 
states (as in the cross term in Eq. (0)): (r 1; r 2 \ip(jimiVij 2 m 2 V2)) a ,i3- The two components of \ip(jimiVij 2 m 2 v 2 )) a ^ 
shown in Eq. (|^) are degenerate in energy. Hence they can interfere with each other, providing interference effects 
that allow for control. 
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To simplify matters we focus on the role of the entangled states |V , C?'i m i' (; :ij2m 2 'i; 2 ))±. Consideration of three special 
cases is in order: 

(a) ji = j%, vi = v 2l mi jt m 2 . Here 

\ip(jim 1 v 1 j 2 m 2 V2))± = — i=(|mi)|m 2 ) ± |m 2 )|mi» ® ® \j2,v 2 )- (13) 

In this case entanglement arises only through molecular rotation motion along a space fixed axis. Case (a) is therefore 
totally analogous to two entangled spins [see Eq. (|l|)]. Indeed, quantum effects arising from such entanglement can 
also be understood in terms of polarization phenomena, in accord with general polarization theory E^ j . 

(b) mi = m 2 ,ji = j 2 , v\ ^ v 2 . In this case we have 

\*{j(jim 1 v 1 j 2 m 2 v 2 ))± = — ;=(|t?i)jV2) ± l^2> l^i)) ® ji) ® |m 2 , j 2 ). (14) 

Here only vibrational motion is entangled; rotational motion is completely separable. Further, unlike case (a), the 
resulting quantum entanglement is entirely due to internal excitation. 

(c) mi = m 2 , vi = v 2 , ji ^ 32- Thus, 

\ip(hmivij 2 m 2 v 2 ))± = -^=(|ji)|j 2 ) ± |j 2 )|ji)) <8> \mi,vi) <g> |m 2 ,u 2 ). (15) 



This case is somewhat complicated. Equation (g_5|) implies that, while vibrational motion is not entangled, rotational 
degrees of freedom arc partially entangled, insofar as the motion along a space fixed axis is still separable. Note that 
in this case quantum entanglement is also due to internal excitation, and that it is intrinsically different from case (a) 
and unrelated to any description based on polarization theory. 

Using Eqs. © and (g), one obtains the following scattering amplitude /±(R) from the initial entangled state 
(k, z, |V , (jimi^ij2m 2 i' 2 ))±) to the final state k', R, j[, m! x , v[,j' 2 , m' 2 , v' 2 \ 

/±(R) = ^f(k''Rj'irn' 1 v' 1 j 2 rn' 2 v' 2 <- kij 1 m 1 v 1 j 2 m 2 v 2 ) 

+^f{k'&j[rn' 1 v' 1 j 2 m' 2 v 2 <- k(-i)j 2 m 2 v 2 j 1 m 1 v 1 ) 

±^f(k / &j[m / 1 v / 1 j' 2 m' 2 v 2 <- fczj 2 m 2 v 2 jimiVi) 

±^/(fe'R-jiwi«ii 2 m 2 ?; 2 <- fc(-z)jiTOiwij 2 m 2 u 2 ). (16) 

Integrating |/±(R)| 2 over R gives |/±| 2 , which is proportional to the cross section for the 3\rn\v' x j 2 w! 2 v' 2 «— 
j\ra\V\j 2 m 2 v 2 transition, with the initial state given by the plus or minus combination. The square of the am- 
plitude for the same transition integrated over R, but from the collision where the initial state is AB(jimit>i) + 
AB(j 2 m 2 f 2 ) [properly symmetrized] is given by |/| 2 . Then, using Eqs. (10) and ( |T6| ) we have that 

J/ + l a + l/-l a . (17) 

By comparing |/+| 2 , |/-| 2 and |/| 2 we can ascertain the effect of entangling rovibrational states, and ascertain the 
degree of control. 

Given Eqs. (|TTJ) , @ and ©, the /±(R) can be written as: 

/±(Il) = T2lW £ £ E V^ l '- i '^' ( R)C^, 2m , 2 ^- mi2 ^ m , 

v v JMm' (j'i2mi2 l'j' 12 m' 12 

x{[l ± (-1)'][C£^^ 

±^j?mr JM (ii<i 2 V 2 il 2 r|i 2U2 Ji«iii 2 0]}. (18) 

Equation (|l^) suggests that there are two important quantum effects in AB + AB scattering associated with 
using entangled molecular states as initial states: (1) Due to the factor [1 ± (—1)'], the permutation sym- 
metry induced by molecular entanglement imposes a parity restriction on the incoming partial waves. That 
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is, for the case of \i/j(jimiVij2m2V2))+ (|V ; 0i m i u ij2^2^2))-), contributions from odd (even) partial waves 
are completely suppressed, whereas contributions from even (odd) partial waves are enhanced. For exam- 
ple, for the case of ^(jiroi^ij^r^i^))-, although the two para-H 2 molecules are spinless bosons, they avoid 
the / — state due to quantum entanglement of rovibrational motion. (2) As indicated by the factor 
[^™ m2 ^ M 0>i^ in Eq. @ , there is quantum 

interference between transitions (iiVij2^2il20 — * (ji v 'ij'2 v 2.fi2^) an d (ji^jxVxjvil) — * {jWxjWij'iaf'') • Whether this 
interference is constructive or destructive depends on the form of the entangled rovibrational state. This effect is signif- 
icant when the magnitude of T JM (j[v' 1 j' 2 v 2 j' 12 l'\jiVij2V2ji2l) is comparable to that of T JM (j' 1 v' 1 j' 2 v 2 j' 12 V\j2V2jiViji2l). 
Further, based on Eq. (|l8|), the scattering amplitude / Qj ^(R) for an arbitrary incoming entangled state 
(fc,z, \ip(jimiVij2m2V2))u,/3) can then be expressed as 

~ cos(qQ + sin(a) exp(^) ^ cos(a) - sin(q) exp(z/3) - 

j a ,p(_ti) -j= j+(ti) h f-W- (J-yJ 

As mentioned above, an alternative but equivalent way to symmetrize the scattering amplitude is to symmetrize 
the outgoing asymptotic state. That is, 

f(k'Ilj' 1 m' 1 v' 1 j' 2 'm' 2 v 2 <— kijimiVij2m2V2) 
= f{k''k,j' l m' l v' l 32m 2 v 2 <- kij 1 m 1 v 1 j2m2V 2 ) 

+f(k'(— R)j2 J Wi^i <~ fczjimivij 2 m 2 v 2 )- (20) 

Based on this procedure, we have 

/±(R) = -^f(k''R.j' 1 rn' 1 v' 1 j' 2 m,' 2 v 2 <- kzj 1 miV 1 j 2 m2V2) 

+ ^f(k\- 1 R.)j'2m' 2 v 2 j' 1 m / 1 v / 1 *- kzj 1 m 1 v 1 j 2 m2V2) 
±^f(k'ilj[m[v' 1 j! i rn2v' 2 <- kzj 2 m2V 2 jim 1 v 1 ) 

±^f(k\-'R)j2m 2 v 2 j' 1 m / 1 v[ <- kzj 2 m 2 V2jim 1 v 1 ). (21) 



In terms of the partial wave expansion, the scattering amplitude /±(R) in Eq. (|2l|) can be further expressed as 

1 i\fn 
V^Vkk 7 



f±&> = E E E V2TTT^'c^ %2m ,c^ 2mi 



JMm' lji? m l2 l'j' 12 m 'i2 
iji2iii2 r / ^iii2 m i2 ym' f-n\rpJM i 



+ ^^^„^" i '(-R)r JM (j2«2Jl«ijl2^|jl«lj2«2jl20] 
±^^^»J^^ m ^™'(R) rJM (^^J2-2^2nj2 V2 Jl« 1 J 1 20 

±^j?w7rVR)r JM (i24il<il2^1j2«2i^iii 2 0]}- (22) 



Interestingly, the equivalence of Eq. ( |22| ) with Eq. ( p.8|) is far from obvious [|19|. As a result, we can no longer readily 
identify the parity restriction condition imposed on the incoming partial waves due to molecular entanglement. In 
this sense, in the theoretical considerations, we prefer Eq. (|lj) to Eq. However, Eq. ( p2] ) is very useful in order 

to confirm the consistency of our results. For all computational examples presented below we have verified that Eqs. 
(U) and (||) give the same results. 

IV. CALCULATIONS AND RESULTS 

In this section we focus on computations for case (c) described above. Specifically, we consider rotational energy 
transfer in rotor-rotor scattering, a problem that has been the subject of considerable interest Jll], 0. We do not treat 
case (a) since it is an example of well-known polarization theory. Case (b), on the other hand, requires computations 
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beyond current capabilities. However, to provide some insight into this case we do present results for one case (b) 
model. That is, we treat a lower dimensional case of vibrational excitation with frozen rotational motion. 

MOLSCAT [^0) is a useful tool for this study. In particular, it provides both cross sections and S'-matrix elements 
necessary for the study of phase control. MOLSCAT also introduces a control parameter to account for whether or not 
the scattering particles are identical. Thus, we can first assume that the two scattering molecules are distinguishable 
and compute all the matrix elements T JM (j' 1 v' 1 j 2 v' 2 j' 12 l'\jiVij2V2ji2l) and T JM (j\ v\ jW?J:?J'' \j2v2j1v1j12l) , and then 
use Eq. (IJ) to calculate the cross sections. The Zarur-Rabitz interaction potential 21 for para-H2 + para-H2 is used 
throughout; it certainly suffices for the demonstration of novel quantum effects in bimolecular scattering. However, 
newer global potentials are available for H2 + H 2 scattering [^2| that can be used for detailed comparisons with future 
experimental studies. 



A. Scattering at E k = 400, 40, 4 cm" 1 

To demonstrate the role of entanglement we consider scattering for different incoming entangled molecular states 
at collision energies Ek = 400, 40, and 4 cm" 1 . In particular, Ek — 400 cm" 1 (« 570° K) is still below the vibration 
excitation threshold, and Ek = 4 cm" 1 (« 5.7° K) is very close to the lowest collision temperature currently achievable 
in molecular crossed-beam experiments. 

We consider the following entangled molecular states in our computations: 

\ipf lh ) = -^(lii>lia) ± b'2>|ii» ® |mi = 0,«i = 0) ® |m 2 = 0,v 2 = 0). (23) 

Subsequent computational investigation |2^] of other a, f3 values in Eq. (||) confirmed that these two states give the 
control extremes. 

Since we consider rotational excitation only, all vibrational quantum numbers are set to be zero, i.e., Vi = V2 = 
v 'i — v 2 — 0- In a typical experiment one might well prepare the AB diatomic superposition state by laser excitation 
of \ji, mi, v\) to \j2,m2,V2)- Assuming linearly polarized light, selection rules ensure that m\ — 1112- However, as 
long as mi = m 2 (so that we can not distinguish between the two scattering molecules by measuring the projection 
of their angular momentum), the results for all m are essentially the same. Hence, here we consider mi = m 2 = 0. 

Figure |l| shows the 9 dependence of the elastic differential cross section a e (9) summed over the final-state quantum 
numbers m'i and m' 2 and integrated over the azimuthal angle <p (R- = (0><£l)j f° r the two incoming entangled states 
\ijjjij 2 ) with ji — 2,j 2 = 0. In particular, Ek equals 400 cm" 1 in Fig. [lp,, 40 cm" 1 in Fig. [l]b, and 4 cm" 1 in 
Fig. [j]c. For all three cases there is a clear difference in a e (0) between the |V^j 2 ) an d iV^j,) cases . Indeed, for 
particular scattering angles, e.g., for 9 close to 7r/2, molecular entanglement induces huge differences. Also evident is 
that the number of minima of a e {9) is always even (odd) for \ipf ■ ) (lip J a ))> a manifestation of the parity restriction 
of the incoming partial waves resulting from the permutation symmetry of the entangled molecular states. Further, 
comparing Fig. [j]c with Fig. |l|a and Fig. it is clear that lower collision energy induces slower oscillations in u e {9). 
This is expected due to the behavior of Y™ (R) in Eq. ( |l8| ) . 

Consider then the total cross section cr e for elastic scattering, denoted a^r for the |V'Jj 3 ) cases. Integrating a e (9) 
over 9 for the lowest energy case (Fig. lc) gives cr+ = 511 A 2 and a~ = 346 A 2 . This is to be compared to the total 
elastic cross section for scattering of AB(j = Q,m = 0) + AB (j = 2, m = 0), proportional to |/| 2 summed over m'i 
and m' 2 , denoted o e and equal to 428 A 2 . If we define the percentage of control in this case as d c = |100(cr+ — a~)/a e \ 
then d c = 39%. For the other two higher energy cases (see Fig. |l|a and Fig. |l|b) the cross sections are found to be 
less sensitive to the quantum entanglement in the incoming state, with d c — 5% and 13%, respectively. Note also 
that the probability of the direct transition j% — 2, J2 = — > j[ = 2, j' 2 = is much larger than that of the exchange 
transition ji — 0, j 2 — 2 — > j[ = 2, j' 2 — 0. Thus, the effect from the quantum interference between the scattering 
amplitudes T JM (20002Z'|20002Z) and T JM (20002Z'|00202Z) [see Eq. @] is negligible. As a result, quantum features 
observed in Fig. |l| are entirely due to the first quantum effect identified in the previous section, i.e., the partial- wave 
parity selection effect. 

In Fig. ||, we show inelastic differential cross sections o~i (9) summed over m'i and m' 2 and integrated over the angle 4>, 
for the incoming entangled states \Tpjij 2 }, with initial-state quantum numbers j% = 4, j 2 = 0, and final-state quantum 
numbers j[ = j' 2 — 2. The collision energy Ek equals 400 cm" 1 in Fig. 40 cm" 1 in Fig. ^b, and 4 cm" 1 in 
Fig. ||c. For these initial-state and final-state channels, the T-matrix elements contributing to the cross sections are 



T JM (2020jy'|40004Z) and T JM (2020jy'|00404Z) [see Eq. @]. A careful examination of these T-matrix elements 
shows that 

T JM (2020j 12 ^'|400040 = (-l) j ™T JM (2Q2Qj' 12 l'\0Q4Q4l). (24) 
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FIG. 1: The 8 dependence of the elastic differential cross section summed over the final state quantum numbers m'i and m'2 
and integrated over 4>. The scattering system is para-H2 + para-H2, and the collision energy is (a) 400 cm -1 (~ 570° K), (b) 
40 cm -1 (~ 57° K), and (c) 4 cm" 1 (~ 5.7° K). Dashed and solid lines are for the incoming free entangled states ) and 

l^j 2 ) [ see Eq - @1> respectively. Here ji = 2, j 2 = 0. 



Hence, the second quantum effect described by Eq. (|i~8|), i.e., quantum interference between T' /M (2020j^'|40004Z) 
and T' /M (2020jy'|00404Z), also becomes important. Indeed, our choice of the final-state channel as j[ = j' 2 is 
motivated by the fact that for j[ ^ j 2 , the second quantum effect is insignificant. As seen in Fig. |^, different incoming 
entangled states cause different oscillatory er^Vpatterns, leading to large differences in <Ji(Q) for fixed scattering angle 
(e.g., for 9 close to or tt in Fig. ||b and Fig. gc, and for 9 around 7r/2 in Fig. ||b). In comparing Figs. ||a, ||b and 
0c, one sees that in the highest collision energy case, quantum effects are least significant. This trend is similar to the 
elastic scattering case, and can also be understood via Eq. (|l8|). In particular, using Eq. (|24|) and the fact that the 
total parity is always conserved, Eq. (^) can be further reduced to 

/±(R) = ~7| J^jj X! V2Z + l(-l)"^ +2 ,, 12 (R)^% n )(- m ; 2 )j; 2 m; 2 Qo%^2m' 1 2m^ 
Jnlj' 12 7n' 12 

xC 4 4 ° 00 [l ± (-l)']T JO (2020 Jl2 (; + 2n)|400Q4/)[l + (-l) l+ H (25) 

where 2n = V - I. Equations (|24|) and (J2|) show that quantum interference between T J0 (2020j 12 r|40004/) and 
T J0 (2020j 12 /'|00404/) can also be interpreted as a selection rule imposed on the quantum number j' 12 - for the incoming 
entangled state \if>4 ) only even incoming partial waves contribute and j[ 2 is restricted to be even, whereas for the 
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incoming entangled state IV^o)' only odd incoming partial waves contribute and j[ 2 is restricted to be odd. This 
is consistent with the fact that the outgoing state is composed of two indistinguishable molecules, with the same 
rotational quantum number (j[ = j' 2 = 2). 





0.0004 Q 2Q 4Q 6Q 8Q 10Q 12Q 14Q 1G[) lgQ 
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FIG. 2: The 9 dependence of the inelastic differential cross section summed over the final state quantum numbers and m' 2 
and integrated over 4>. The scattering system is para-H2+ para-H2, and the collision energy is (a) 400 cm -1 , (b) 40 cm -1 , and 
(c) 4 cm -1 . Dashed and solid lines are for the incoming free entangled states \ipt.j,} and |^^j 2 ) [see Eq. (Ba)], respectively. 
Hereji =4, j 2 =Q,j[=j' 2 = 2. 

The case in Fig. ||c is worthy of further discussion. In this case one sees that 0^(6*) for \^J x j 2 ) (solid curve) is 
systematically smaller than that for |^^j 2 ) (dashed curve), for almost the entire range of the scattering angle. This 
leads to a significant difference in the total inelastic cross section cr- for the \%bf ■ ) cases: &f = 0.057 A 2 and 

= 0.032 A 2 , giving a control range of d c — |100(<7^ — a i )/<Ji)\ =57%, where cr^ is the total inelastic cross section 
for scattering from AB(j = 4, m = 0) + AB (j = 0, m = 0) to the channel j[ = j' 2 = 2. However, d c is less than 2% 
at the higher energies in Fig. |[ Thus, as in elastic scattering, total inelastic cross sections can also be sensitive to 
the quantum phase embedded in initial entangled molecular states at sufficiently low energies. 



B. Ultracold Molecular Collisions: Extreme Quantum Effects 



Given the results for Ek > 5° K, we consider ultracold collisions. This new area of focus, ultracold molecular 
scattering, has already displayed interesting quantum effects [p4[. Further, understanding this area is also crucial in 
realizing molecular Bose-Einstein condensation 



25, M 
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From the control point of view, the advantage of ultracold collisions is that scattering cross sections are often 
dominated by very few partial waves. This feature of ultracold collisions may well allow for new opportunities in 
manipulating molecular collisions ^7|. For example, the fact that entangled incoming molecular states impose parity 
restrictions on the partial waves [see Eq. (fl8|) 1 implies that quantum entanglement should be even more important 
at ultracold temperatures. Consider, for example, a zero collision energy limit case where only s-wave scattering 
(/ = partial wave) contributes to the cross sections p8| . In this case, if s-wave scattering is selected by molecular 
entanglement, the result would be enhanced nonzero cross sections. If s-wave scattering is forbidden due to molecular 
entanglement, then molecular scattering is completely suppressed. That is, in the elastic scattering case, all the 
scattering would be in the forward direction, and in the completely suppressed inelastic scattering case, there would 
be no inelastic scattering in any direction. 

To confirm such theoretical considerations and to determine how cold the collisions need be in order to experimen- 
tally observe such extreme quantum effects, we extend our calculations to the ultracold regime, with Ek ranging from 
0.4 cm -1 to 0.0004 cm -1 . These computations are generally easier than at higher energies because the maximum 
total angular momentum contributing to cross sections is very small (29). 

Figure || shows a e {9) for two energies [0.4 cm -1 and 0.04 cm -1 ], with the incoming asymptotic states given by 
l^i,) (ji — 2, 32 — 0). In both cases a e (6) for |V>t,-„) (dashed line) is fairly uniform suggesting that s-wave dominates 
the elastic scattering. Thus, preventing s-wave scattering via entanglement is expected to cause a dramatic decrease 
in the cross section. Indeed, in Fig. ||b, one sees that elastic scattering is almost totally suppressed for the \ipZj 3 ) 

case. Integrating cr e (9) over 9 gives er+ = 724 A 2 and a~ = 285 A 2 with d c = 87% for the case shown in Fig. |^a, and 
cr+ = 1014 A 2 and a~ = 11 A 2 with d c = 195% for the case shown in Fig. |b. 




FIG. 3: As in Fig. [l] except that elastic collision is ultracold. The collision energy is (a) 0.4 cm 1 , and (b) 0.04 cm I . 

Figure [| shows sample results for ultracold inelastic scattering. As in Fig. ||, IV'jj*,) (ji — 4, j% — 0) is chosen as 
the incoming asymptotic state and we examine transitions to j[ = 2,j' 2 = 2 with collision energies of 0.04 cm -1 and 
0.0004 cm -1 . The differences in <Ji(9) between |V$[j 2 ) an d \'<Pj 1 j 2 ) m Fig. |^a is somewhat similar to that in Fig. ||c, 
where the collision energy is 100 times larger, with a corresponding d c of 22%. Hence, control is significant, but the 
system is insufficiently cold to show extreme quantum effects. By contrast, Fig. ||b shows results for the lowest Ek of 
0.0004 cm -1 . Clearly, in this case, inelastic scattering for | - 0^j 2 ) 1S a l mos t completely suppressed, whereas inelastic 
scattering for [V^^) remains significant. The integrated cross sections give af = 1.63 A 2 and err = 0.02 A 2 , with a 
4=194%. 

Thus, we have demonstrated that quantum entanglement in the incoming asymptotic state can dramatically alter 
ultracold collision cross sections. These effects will also appear in the scattering of the superposition states in Eq. 
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(§). That is, Eq. (|) is well obeyed in all of the cases shown in Figs, [j] to ||. For example, in the elastic scattering 
case (ji = 2, j 2 = 0, vi = v 2 = 0, mi = m 2 = 0, v[ = v' 2 = 0, j[ = 2, j' 2 = 0) the scattering amplitude for the 
ji = 2,j2 = 2 — > j( = 2, 72 = transition due to the first satellite state is orders of magnitude smaller than that for 
the ji — 2,j2 = — * j[ = 2,j 2 = transition associated with the entangled state term at the Ek considered. The 
probability of the transition j\ = 0,j2 = — > j[ = 2,j 2 = due to the second satellite state is also negligible (indeed, 
it is exactly zero for Ek < 364.8 cm -1 ). Similarly, for the inelastic scattering case (ji = 4, j 2 = 0, v\ = v% = 0, 
mi = m 2 = 0, v[ — v 2 = 0, j[ = 2, j' 2 = 2), the scattering amplitude due to the first satellite, the de-excitation 
ji = 4, j2 = 4 = 2,j 2 = 2, is negligible due to the large rotational energy mismatch between the initial and 
final channels. In addition, the second satellite state does not contribute to the product channel (j[ = 2, j 2 — 2) for 
E k < 729.6 cm -1 . 
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FIG. 4: As in Fig. g except that inelastic collision is ultracold. The collision energy is (a) 0.04 cm , and (b) 0.0004 cm . 
Thus, for example if, in Eq. (|J), a\ — a 2 = 7t/4, then Eq. (||) reduces to 

<?{0Md P )) w \<r(9,\i>+ h )), P = 2mr; 

<r(0,\1>dp)) « \°{0Mj lh )), /9=(2n+l)7r. (26) 

In this case the role of (3 as a control parameter, evident in Eq. ( |26j ) as a method for switching between the plus and 
minus states, is worthy of note. For example, the elastic scattering results in Fig. show that by varying (3 one 
can enhance, or almost completely suppress, the yield of the (j[ — 2,j 2 = 0) channel in directions other than that of 
the incident momentum. Likewise, the inelastic scattering results in Fig. ^b suggest that by manipulating (3 we can 
enhance or almost shut off the (j[ — 2, j' 2 — 2) channel. 

C. Vibrational Relaxation 

The above discussion provides results on rotational energy transfer. The formalism, however, applies equally 
well to other types of scattering. Unfortunately, the current state-of-the-computational-art prevents a full quantum 
calculation on reactive scattering, or even on diatom-diatom rovibrational energy transfer. That is, the exact numerical 
treatment of AB + AB in three dimensions is still not feasible and approximations such as the sudden approximation 
may not be appropriate for low energy collisions of the type that we are examining since it assumes that molecular 
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rotation is slow compared with that of vibration and translation. Further, although a semiclassical treatment of 
rovibrational energy transfer is available |3l| it ignores the role of the quantum effects in the translational motion. 
Nonetheless, we wish to gain some insight into phase control of vibrational relaxation. For this reason we present 
results on a simplified model of diatom-diatom scattering where the rotational motion is frozen. Specifically, we 
adopt a three-degree-of-freedom model that assumes that both the projectile and target molecules point at a fixed 
direction during the entire scattering process. The scattering problem is then computationally solved by wavepacket 
propagation. To do so we employ a time-dependent approach with real L 2 eigenfunctions with damping p2, 
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5: The inelastic transition probability of the channel (v[ = 1, v' 2 = 1) versus the collision energy in the para- H2 + para- 



H2 scattering, for three incoming entangled vibrational states described by Eq. (|27|). a = tt/A, /3 = (top curve), n/2 (middle 
curve) and ir (bottom curve). The internuclear axis of each para-Hb molecule is assumed to be parallel to the incident velocity 
during the scattering process. 



Figure |5| displays the transition probability to the product (v[ = l,v' 2 = 1) channel for three initial entangled 
vibrational states characterized by 

= cos(a)|«i = 0) ® \v 2 = 2) + exp(i0) sm{a)\vi = 2) <g> \v 2 = 0) (27) 

with a = 7r/4, (3 = 0, tt/2, and ir, respectively. As seen in Fig. ||, the quantum phase j3 embedded in the initial en- 
tangled vibrational states strongly affects vibration- vibration relaxation, with the cross section changing considerably 
as a function of the initial entangled state. 



V. DISCUSSION AND SUMMARY 



The results on para-H2 + para-H2 scattering make clear that entanglement can play a significant role in the 
control of nonreactive collisional processes. To investigate this phenomenon experimentally one could attempt to 
prepare the entangled molecular state directly. However, although quantum entanglement of atomic systems has 
been experimentally realized |34|, |3^] , preparing entangled molecular states would require a considerable extension of 
technology. By contrast, contributions from entangled states appear naturally in the collision of AB + AB systems 
where the initial states are prepared as superpositions of the form given in Eq. (|^). Preparing these states would 
entail excitation of the \j\,mx,Vx) to produce the superposition with lj'2, wi2, i>2)- If the transition is dipole allowed 
then direct excitation of the lower state is possible. Alternatively, stimulated Raman adiabatic passage (STIRAP) 
|56| or the tripod-STIRAP scheme |37| may provide a useful choice for preparing the superposition states for both 
the projectile and the target molecules. In the latter approach an additional laser couples the intermediate level 
(through which the initial and final state are radiatively connected via the pump and Stokes-laser) with another 
unpopulated state. Depending on the time overlap of the additional laser with that of the pump and Stokes-laser, 
any coherent superposition state of the initial and final state can be created. Thus, by introducing different coherence 
characteristics from laser fields into the scattering system, we can control (3 and thus select the form of the entangled 
state component \ipdp) m Eq. (^), giving rise to phase control over differential and total cross sections. 

Our studies have been restricted to nonreactive scattering of AB + AB, with particular numerical application to 
para-H2 + para-H2. As such, it remains to discuss the potential for applications to inelastic scattering of other 
systems and to comment on control of reactive scattering of identical particles. 

General inelastic AB + AB scattering: Extensions to other inelastic scattering cases fall into two categories. The 
first are systems which, like para-H2 + para-H2, have zero total nuclear spin, zero total electronic orbital angular 
momentum, and zero total electronic spin. The theory of these cases is described above and is directly applicable 
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to inelastic scattering of other zero spin cases such as CO + CO, a system whose vibrational energy transfer is of 
interest in laser physics. However, the choice of zero spin merely simplifies theoretical considerations and it is not 
essential to control. 

The second category, molecules with nonzero nuclear spin or electronic spin, does require an extension of the theory 
provided above. For such molecules, one of the two quantum effects due to entangled molecular states, i.e., the 
parity selection of the incoming partial waves, may be substantially reduced when the nuclear or electronic spin is 
unpolarized. To see this note that the spin degree of freedom can introduce additional permutation symmetries into 
the system. Thus, entangled rovibrational states like \ip(ji v i m ~tj'2 rn 2 v 2)}± do not necessarily select the partial waves. 
For instance, for two unpolarized identical molecules with total integer nuclear spin I/v, the permutation symmetry of 
all other degrees of freedom is +1 with the probability (Jjv + l)/(2Ijy + 1) and —1 with the probability (In)/(2In + 1) 
. When J/v is large, two permutation symmetries of all other degrees of freedom are equally allowed. As a result, the 
permutation symmetry induced by entanglement of molecular rovibrational motion will have a negligible effect on the 
differential or total cross sections. However, as suggested in Ref. fHj], one quantum effect still survives, i.e., quantum 
interference between the transition amplitudes T JM (j^v^j^j'^l uivij2V2ji2l) and T JM (j' 1 v' 1 j 2 v 2 ji2l'\j2V2jiViji2l)- 
Establishing the magnitude of this effect alone will require further study. 

Finally, note that there are considerable differences between the inelastic scattering case studied here and the case of 
reactive AB + AB scattering (to form A2B + B or A2 + B2), the subject of future research. A full treatment of reactive 
scattering, which does require considerable extensions of current computational technique, is under consideration. 

In summary, we have examined the nature of the interference, and hence control, in nonreactive AB + AB scattering 
when the initial states are in a quantum superposition state. In doing so we have exposed the relationship of the 
interference term to entangled molecular rovibrational states as the incoming asymptotic state. Intriguing quantum 
effects resulting from quantum entanglement are revealed, and are shown to provide a novel means of controlling both 
nonreactive differential and total cross sections in identical diatom-diatom scattering. 
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